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Abstract 

A model for the coexistence of p-w&ve superconductivity (SC) and ferroelectricity (FE) is pre- 
sented. The Hamiltonian of SC sector and FE sector can be diagonalized by using the so(5) and 
h{4) algebraic coherent states respectively. We assume a minimal symmetry-allow coupling and 
simplify the total Hamiltonian through a double mean-field approximation (DMFA). A variational 
coherent-state (VCS) trial wave-function is applied for the ground state. It is found that the fer- 
roelectricity gives rise to the magnetic field effect of p-wave superconductivity. 
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I. INTRODUCTION 



The coexistence of superconductivity (SC) and ferroelectricity (FE) and competition play 
a important role in motivating the work of Bednorz and Muller [l| on the high tempera- 
ture cuprates superconductors. Even earlier, work on the "old" superconductors of the (3-W 
structure like V3S1, Nb 3 S n etc., where T c ~ 23K and a martensitic phase transition occurs in 
the same temperature range, gave rise to investigations on the possibility of a "ferroelectric 
metal" or a "polar metal" and thus to the study of SC-FE coexistence 0, Q- S- 0- El Many 
theoretical papers have already studied microscopic models for the effect of lattice insta- 
bility on superconductivity in the sodium tungsten bronze systems 0, |], J These papers 
have illuminated many aspects of the interplay between the structure deformations, such as 
rotation of underlying octahedaral units and coupling with electron pairs. Study of SC-FE 
coexistence problem can be relevant to recent work by Weger and collaborators jloL , on 
the mechanism of high temperature superconductivity in the cuprates. In that work the 
presence of a nearby FE instability close to the SC transition is related to the anomalously 
large ionic dielectric coefficient in the cuprates, which reduces the electron-electron repul- 
sion and then can lead to an enhanced net electron-electron attraction, producing higher 
T c . Moreover the recent so(5) and su(A) models for multi-critical superconductor anti- 
ferromagnetic behavior in the high temperature superconductors have also been studied by 
Zhang et.al. jl^ . ^5|. On the other hand, the recent experimental discovery of the 

coexistence ferromagnestism (FM) and SC in VGe^ an d subsequently in ZrZn 2 |3| 

and URhGe ^J], have shown clearly that the spin-triplet states of p- wave SC are realized 
in the nature. This naturally causes our interest in the relationship between FE and p-wave 
SC. This paper is organized as follows: In Sec. II we will construct a so(5) <g> h(4) algebraic 
structure general model for the coexistence of p-wave SC and FE. In Sec. Ill and IV we will 
diagonalize the Hamiltonian of p-wave SC and FE according to so(5) and h(4) algebraic co- 
herent state methods respectively, and get their energies, the eigenstates and the expection 
value of the order parameters in the ground coherent state. In Sec. V the total Hamiltonian 
including the biquadratic interaction is simplified to the bilinear Hamiltonian H^mfa by 
using a double- mean- field approximation. We can find the variational solution of H DMF a 
by forming a trial eigenstate analogous to the product coherent state. The energy spectrum 
and the eigenstates will be obtained based on the variational solution. 
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II. OUR MODEL 



Here we will consider a general model for the coexistence between p-wave superconduc- 
tivity and ferroelectricity. According to a known result there is so(5) structure in p-wave 



superconductivity [2JJ, |2l| that is formed by two su(2) not commuting with each other, 
where one describes the attractive BCS interaction and the other the usual spin opera- 
tors, as well as other four generators associated with the transitions. But ferroelectricity is 



formed by Heisenberg algebra h(A) (22J, |23j, |24j, |25(. Motivated by Ref. |26(, we write the 



total Hamiltonian in three parts as follows: 

H = H S c + Hfe + H INT , (1) 

where the first term Hsc is the BW type of p-wave superconductivity mean-field reduced 
Hamiltonian, namely 

H SC = Y, h ^ (2) 

k 

with the Hamiltonian at each given momentum k given as 

K = e k ^ k) + [A E (k)E { ^ + A [/ (k)[/| k) + A v (k)vl k) + H.c.}. (3) 
Here A E (k) = \ £ Vyx < E [K) >, A [/ (k) = i£^kk' < > and A y (k) = \ E Vkk' < 

k' k' k' 

are opposite and equal spin pairing ( "gap" ) energy respectively, and j)- wave attrac- 
tion pair interaction potential = —3Vi(k, k')k ■ k'. 

The second term H FE is the displaced oscillator Hamiltonian of the displacive ferroelectric 
soft mode phonon |2T]. It's taken as 

H FE = u TO {b ] b +\) + li^ + b), (4) 

where uto is the frequency of the soft TO mode, 71 is taken as a positive constant, and e is 
the magnitude of the electric field E. 

The third term Hint is the interaction term. According to Ginzburg-Landau (GL) theory 
every term in the free energy shall be a scalar invariant under the relevant symmetry 
group, here the free energy density of system will be the form of 5F = a\P\ 2 + ||P| 4 + 
a|A| 2 + /5|A| 4 -|-fi;|A| 2 |P| 2 . In such a GL theory, the lowest-order, generic coupling between 
superconducting A and the ferroelectric P will be a biquadratic term of the form proportional 
to |A| 2 P 2 , that is, the last term in SF. This lowest-order term satisfies gauge and parity 
symmetry requirements. From the correspondences |A| 2 ~ (E + E- + U + U- + V + V- + H.c.) 
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and P 2 ~ (6 + tf) 2 we can immediately translate this coupling term into the interaction term 
of our model as follows: 



H INT = ^ 72k (4 k) Ei k) + U^U { y + vfVi k) + Hc)(tf + b) 2 , 



(5) 



where 72k is the initial pair-TO mode coupling coefficient. 

In our model, the set {E±\ F^ k \ U± \ V±^} form so(5) algebra, its generators 

are expressed as 



E + ] - 72( a -kT a k| + a -k| a kt)' - 75 ( a k| a kt + a!_ k | a -kT)> 



U+ ^ — al k t a kt' 



^+ k) = a -k|°ki ? v -' = Ok|a_ ki , 

^3 = |( a kT a kT + a -kt a -kT + a k|°k| + fl tk| a -k| — 2), 



77 (k) - n n 

u_ — ak|a_kt) 

Kk) 



(4| a k| + al k | a -k| - 4T a kt - al k T a -kT)> 
and satisfy the following commutation relations: 



[4 k) 


,^ k) ] 


= ±Ff, 




= TE?\ 




= ±4 k) 


[4 k) 




= ±4 k) , 


[4 k) ,4 k) ] 


= TVi k \ 


[4 k) ,Ff] 






,4 k) ] 


= ±4 k) , 




= ±f£\ 




= ±U^ 






= TU?\ 




= ±vf\ 




= ±vi k) 


[4 k) 


X k) ] 


- P( k ) 

— ^3 > 


[uP,uM] 


= (Ei k) - 






[F (k) 


)F (k)] 


- ^3 > 


[vi k) M k) ] 


= (4 k) + 







(6) 



(7) 



Note that Hsc is linearized with respect to so(5) algebraic generators for a given momentum 
k, Hfe is written in terms of h(4) generators, and Hint is expressed as a product of the 
bilinear form of so (5) and quadratic form of h(A). Therefore the total Hamiltonian is a 
so{h) <S> h(4) direct product algebraic structure in each given momentum k. 



III. THE so(5) STRUCTURE OF p-WAVE SUPERCONDUCTIVITY MODEL 

From Hamiltonian Hsc we can known that is written in terms of so(5) generators for 
a given k. The dynamical symmetry or spectrum generating algebra for each k is so(5)k, so 
the spectrum generating algebra of Hsc is ®k<so(5)k- 
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The eigenstates of H S c are expressed by a direct product of so(5) k coherent states <8> k |£ k )- 
Therefore the eigenstates |£) can be written as 



where 



IO = ®k|&> = ®kW^k)b,g>, (8) 

W(f k ) = exp{^ k ( v / 2cos^ kJ Bi k) - sin fl k e^ k f/[ k) 
+ sin^ k e-^ k Vi k) ) -H.C.}, 

with the real coherent parameters £ k , (#k, 0k) are angles in spin space for a given momentum 
k. \p,q) is the mutual eigenstates of the Cartan subalgebra {E^\ F^} of so(5), and their 
eigenvalues are p and g respectively. By tedious calculations we can immediately diagonalize 
the Hamiltonian h k as 



W\^)KW{^) = y^ + A^k)^, (9) 
where A 2 (k) = 2(|A E (k)| 2 + |A[/(k)| 2 + |A y (k)| 2 ). Also we obtain the gap equation: 



A B (k) 
At/(k) 

Ay(k) 



kk' 



A(k') 



^§ cos 6» k / 

^sin^e^' 

^±2) S in0 k ,e-^' 



(10) 



From the above gap equation we can obtain the following results: 

1. If both p and g are zero, the gap A = and the p-wave SC lies in disorder state. Its 
eigenstate is \£ dis >= ® k W(£ k )|p = 0, q = >. 

2. If either p or g are not zero, the gap A ^ and the p-wave SC lies in superconducting 
state, but when p — — 1 and q = the p-wave SC lies in the ground state. Its ground state 
is \Cbcs >= <S>kW(£k)b = — 1, q = > where |p = —1, g = > is the vacuum state. 



IV. A h{4) STRUCTURE OF DISPLACIVE FERROELECTRIC SOFT MODE 

The Hamiltonian H FE can be transformed by the displaces oscillator Bose operator 

f/(Co)=exp[Co(& t -&)], (11) 



where the coherent parameter ( is taken as real, so we obtain 



(lie? 



(12) 



Note that U j H FE U is shifted to a new minimum, but it retains the same excitation frequency 
ojto as the original oscillator. The eigenstates and eigenvalues of Hfe can be given as 



|Co >= U(( )\n >= exp[Co(6 t - b)]\n >, 



(13) 



and 



ijJTo{n + 



(7i£) 



(14) 



2 uto 

where \n > is a number eigenstate of the phonon number operator 7V& = tfb and the state 
| Co > is a Glauber coherent state for the FE oscillator 29l l30|. 

The order operator for the FE polarization is the coordinate operator Q or (b + ¥), 
therefore the order parameter 



Vfe 



<Co|& t + &|Co>=2Co 
< n\b ] + b\n >= 0. 



When rjFE = it shows that the phonon is free, but tjfe — 2Co h exhibits that the phonon 
is polarized spontaneously. 



V. THE VARIATIONAL COHERENT STATE EIGENSTATES OF OUR HAMIL- 
TONIAN UNDER DOUBLE-MEAN- APPROXIMATION 

The total Hamiltonian is turned into bilinear forms by double mean field approximation 
procedure which reduces biquadratic operators such as A 2 B 2 to the following form 

A 2 B 2 w (2 A < A > - < A > 2 ){2B <B> - <B > 2 ), (15) 

based on the assumption (A— < A >) 2 m and {B— < B >) 2 m 0. After making the 
double-mean-approximation and isolating a single mode k, we have the effective Hamiltonian 
at each momentum mode k in the DMFA: 

H DMFA = eE 3 + (A' E E+ + A' V U + + A' V V + + H.c.) + u TO (tfb + -) + ^(6+ + b) 

+(rf E + + Y U 2 U + + Y V 2 V + + H.c.)(tf + b) + T 3 . (16) 
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Therefore the double mean field approximation (DMFA) yields a bilinear effective interaction 
term, and it renormalizes the coefficients A' and ^2,3 in Homfa as follows: 



A E = A E - 72 < 6 f + b > 2 < E_ >, 
A^ = A v -72 < 6 f + b > 2 < U- >, < 
Ay = A y - 72 < 6 f + b > 2 < V_ >, 



rf = 2 72 <tf + b>< E. >, 

If = 2 72 < 6 f + b >< U- >, 

r 2 r = 2 72 <6 t + 6><\/_ >, (17) 



Tx = 7l e - 2 72 < 6 f + b > (< >< > + < [/+ >< [/_ > + < V+ >< VI >), 
r 3 = 72 (< >< > + < U+ >< U- > + < V+ >< VI >) < 6 f + b > 2 . 

This total Hamiltonian includes the pure p-wave SC and the pure FE prototype systems 
and their coupling via the soft-mode oscillator coupled to opposite and equal-spin pairing 
Hamiltonian. Noting that the initial Hamiltonian is the enveloping algebra of so(5) ® h{A) 
because of the biquadratic interaction terms, while H DMFM is an element in the direct 
product algebra so(5) <8> h(A). 

When 72 — > 0, we recover the sum of two separate sectors: for p-wave SC and FE. In 
order to obtain the ground state eigenstates and eigenvalue of our model HdmfAi we will 
make use of the variational principle. Introducing an analogous trial variational coherent 
state (VCS) which is the product of two coherent-like states and is denoted \ip v >: 

\<p v >= |f > |C >= V ± \p, q>V 2 \n >, (18) 

where Vi(f) = exp[f (V2cos9E + -sm9e^U + + sm9e-^V + )-H.c] and V 2 (C) = exp[C(& t -6)]. 
Here the real parameters f and ( are variational unknowns. The kets \p,q > and \n > are 
the same as before. Now we define the energy in state \ip u > as the diagonal value of H^mfa 
in the variational coherent state \ip v >: 

= pe cos 2f + [--^= A E cos + V -^- A^e"'* sin 9 - V -^- A y e^ sin 9 + H.c] sin (2f) 
+2C[--^rf cos 9 + P^-rUe-i* S in 9 - ^t^e** sin + H.c] sin (2f) 

+c^ TO (n + c 2 + 1) + 2riC + r 3 

= pe cos (20 + pA sin (2f ) + cu TO (n + C 2 + ^ 

+672CV + g 2 sin 2 0) sin 2 (2f) + 2 7l eC. (19) 
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Here the coefficients A' and 1^2,3 in the energy E Ptqtn (£, 9, (j), () are given as 



A E = A E - 4 72 C 2 sin (20 (-^ cos theta) = -j- cosfl[A - 4 72 C 2 psin (20], 
A' v = A c/ -4 72 C 2 sin (20(2^ sin 9e**) = \ sin 9e^[A - 4 l2 ( 2 (p - q) sin (20], 
A' T/ = A v -4^ 2 sinf2nf- 2 ± 2 sin# e -^) = -\ sin 0e~^[A - 4 72 C 2 (p + q) sin (20], 



1 2 



1 2 



r - A y -4 72 C 2 sin(20(- £ ?sini 
= 4 72 Csin(20 



1 2 



p-q 



sin I 



-2±2sini 



I\ = 7l e - 2 72 C(p 2 + q 2 sin 2 0) sin 2 (20, 
T 3 = 2 72 C 2 (p 2 + q 2 sin 2 9) sin 2 (20- 

We determine 9, 0, £, £ from the following equations: 



(20) 



8E/d9 = 


o. 


dE/d(f> = 


0, 


dE/di = 


0. 


dE/dC = 


0. 



(21) 



Eqs. (|2~T|) have two solutions: 
sin 9 = 0, 

tan2£ = [A + 12 72 C 2 psin (2£)]/e, 

. e = ~[7ie + 6 72 CP 2 sin 2 (20]M-o- 
This result exhibits that FE and SC can coexist in the opposite spin pairing state of 

p-wave SC but the equal spin pairing state of p-wave SC disappeared. Here the energy is 

E = pe cos (20 + pA sin (20 + uj TO (n + C 2 + ~) + 6 72 C V sin 2 (20 + 2 7l eC- 

cos 9 = 0, 

tan 2^ = [A + 12 72 C 2 (p + f ) sin (20]/e, 

. £ = -[Tie + 6 72 C(p 2 + q 2 ) sin 2 (2£)]/uj to . 
This result exhibits that FE and p-wave SC may coexist in the equal spin pairing state 

(ABM state) but the opposite spin pairing state of p- wave SC disappeared. Here the energy 



is E = pe cos (20 + pA sin (20 + oj TO {n + C 2 



6 72 C V + q 2 ) sin 2 (20 + 2 7l eC 



Therefore from above two cases we get that the spin of p-wave SC is oriented because 
of the existence of FE. It indicates that the FE effect in the excite state of p-wave SC is 
equivalent to magnetic field. 
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VI. CONCLUSION 



In this paper we have constructed a general model for the coexistence of p-wave su- 
perconductivity and ferroelectricity. The Hamiltonian of p-wave SC can been diagonalized 
based on so(5) spectrum-generating algebra structure and has shown that the eigenstate is 
related to so(5) coherent state. Also the Hamiltonian of FE is diagonalized by using h{A) 
algebraic coherent state. The total Hamiltonian under the double mean-field approximation 
can been solved by making use of a variational coherent-state (VCS) procedure by virtue of 
the so(5) <S>h(4) algebraic structure in each given momentum k. This leads to the conclusion 
that the ferroelectricity gives rise to the magnetic field effect of p-wave superconductivity. 
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